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Abstract. We present UV rest frame spectra of 3 power- 
ful narrow line radio galaxies and the hyperluminous 
type 2 active galaxy SMM02399-0136, all at high redshift 
>2). We find high velocities (FWHM > 1000 km s" 1 ) 
in the extended gas of all objects. A natural explanation is 
the interaction between the radio jet and the ambient gas, 
that drives shocks into the gas and accelerates the clouds. 
However, the existence of high velocities in regions where 
such interactions are not taking place implies that other 
processes can play a role. We discuss here several possible 
mechanisms. 



X 



1. Introduction 

The existence of high velocities (FWHM>1000 km s _1 ) 
(McCarthy et al. 1996) in the extended gas (EELR) of 
high redshift (z >2) radio galaxies (HZRG) is in con- 
trast with the more relaxed kinematics observed in the 
majority of low redshift radio galaxies (FWHM < 400 km 
s _1 ) (Tadhunter et al. 1989). The nature of such extreme 
kinematic motions is not well understood. We investigate 
here this issue studying the kinematics of the extended 
gas in a small sample of 4 distant active galaxies (z >2): 
MRC1558-003, MRC2025-218 and MRC2104-242 (radio 
galaxies) and SMM02399-0136 (hyperluminous type 2 ac- 
tive galaxy with very weak radio emission). 

2. Observations and data reduction 

The spectroscopic observations were carried out on the 
nights 1997 July 3-5 and 1998 July 25-27 using the EMMI 
multi-purpose instrument at the NTT (New Technology 
Telescope) in La Silla Observatory (ESO-Chile). The de- 
tector was a Tektronix CCD with 2048x2048 pixels of 
size 24 /im, resulting in a spatial scale of 0.27 arcsec per 
pixel. We used EMMI in RILD spectroscopic mode (Red 
Imaging and Low Disperson Spectroscopy). We used the 
same grism (#3) for all objects. This has a blaze wave- 
length of 4600 A, dispersion 5.9 A/pixel and wavelength 



Send offprint 
lard@iap.fr 



requests to: M.Villar-Martm. email 



range 4000-8300 A. The slit was aligned with the radio axis 
for the three radio galaxies. We positioned the slit along 
the two main optical components of SMM02399-0136 (LI 
and L2, adopting the nomenclature of Ivinson et al. 1998 
[IV98 hereafter]). A log of the spectroscopic observations 
is shown in Table 1. 

Standard data reduction techniques were applied using 
IRAF software (see Villar-Martin et al. 1998 for a more 
detailed description) . 

3. The fitting procedure 

In order to study the kinematics of the gas, we fitted 
the emission lines with a Gaussian profile at every spatial 
position (pixel). Several spatial pixels were added where 
the emission was too faint. We used the Star link packa- 
ge DIP SO for this purpose. The FWHM, flux and central 
wavelength were measured from the Gaussian fitted to the 
line profile. The FWHM was corrected in quadrature for 
instrumental broadening (the instrumental profiles in the 
observed frame are given in Table 1). 

Single Gaussians did not always provide a perfect fit 
and underlying broad wings were sometimes present. This 
is probably due to the presence of several kinematic com- 
ponents and/or absorption of Lya by neutral hydrogen. To 
eliminate uncertainties due to the second mechanism, we 
also present the result of the fit for the second strongest 
emission line CIVA1550 not susceptible of hydrogen ab- 
sorption. 

At the s/n of the data, we are confined to using sin- 
gle Gaussian fits to the lines, a procedure which is the 
same as that followed by Tadhunter et al. (1989). A single 
Gaussian fit will a) lose any information about multiple 
components b) neglect any possible weak broad under- 
lying wings (as is observed in MRC2104-242, see Fig. 1 
left cloud), since the fit will be optimized for the domi- 
nant part of the line. However, the main goal of this paper 
does not require such a precise analysis. The broad wing 
on MRC2104-242 will have little influence on the fit, which 
will be dominated by the strong, narrower component. 
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Table 1. Log of the spectroscopic observing run. The redshift (calculated from the spatially integrated Lya emission) 
is also indicated. 
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MRC2 104-242 


2.492 


3/Jul/97 


4x1800 


1.5 


1.10 


10.4 


22 


MRC2025-218 


2.632 


3/Jul/97 


3x1800 


1.5 


1.15 


10.9 


30 


MRC 1558-003 


2.530 


25-26/Jul/98 


4x1800 


1.5 


1.00 


10.3 


75 


SMM02399-0136 


2.803 


25-26/Jul/98 


5x1800 


1.5 


1.20 


10.3 


88.6 



4. Results 

We present in Figs. 1 to 4 the results of our analysis for the 
4 targets in the sample. The upper panel in each figure is 
the 2-D spectrum of the Lya spectral region with the dis- 
persion in A running vertically. The middle panel presents 
the spatial variation of the FWHM and the bottom panel 
shows the spatial variation of the velocity shift of the Lya 
emission (open circles) across the nebula. The 3 panels in 
each figure have the same spatial scale and are aligned so 
that vertical lines join the same spatial positions. FWHM 
for CIVA1550 is also plotted for comparison (solid trian- 
gles). 

- MRC2104-242 (Fig. 1): Lya shows a bimodal distri- 
bution and is extended over ^12 arcsec along the slit 
aligned with the radio axis. The two blobs lie in be- 
tween the radio lobes (McCarthy et al. 1990). The bi- 
modal distribution is also apparent in CIV. The two 
clumps present high FWHM values (1100 km s _1 and 
900 km s _1 respectively, for the spatially integrated 
spectra) and are shifted by ^500 km s _1 (consistent 
with McCarthy et al. 1990, Koekemoer et al. 1996). 
Kinematic substructure is observed in the two blobs. 
The velocity curve (Fig. 1, bottom panel) is rather flat 
across each blob. CIV is also extended and presents 
high FWHM values (-700 km s _1 , see Fig. 1). 

- MRC2025-218 (Fig. 2): Lya shows a bimodal distribu- 
tion and is extended over ~4.5 arc sec. This structure 
lies between the radio lobes (Pentericci et al. 1998). 
Continuum is detected in our spectra. Large and rather 
constant FWHM values are measured across the two 
components (^1200 km s _1 and 700 km s _1 respec- 
tively). The velocity curve is rather steep across the 
nebula varying smoothly over a range of ~600 km s _1 . 
CIV is also extended and presents similar FWHM as 
Lya. 

- MRC1558-003 (Fig. 3): Lya is extended over >15 arc 
sec. A bright component is detected as well as diffu- 
se, very extended emission which show large velocity 
widths at ~10 arc sec from the main component. The 
optical (Rotggering et al. 1994) and radio astrome- 
try (Rottgcring et al. 1996, Rhee et al. 1996) locate 
the high velocity region several arcsec beyond the ra- 
dio structures. CIV is also extended and presents very 
large FWHM within the high velocity region, consis- 



tent with the Lya measurement. There is no apparent 
pattern in the velocity curve of the ionized gas. 

- SMM02399-0136 (Fig. 4): This object is a hypcrlumi- 
nous active galaxy, gravitationally lensed by a fore- 
ground cluster (IV98). It consists of two main optical 
sources LI and L2. The radio emission is very weak, 
below the detection tresholds of most radio surveys. 
The radio, submm and optical properties are consis- 
tent with a scenario where LI contains an active nu- 
cleus and L2 is an interacting companion. The system 
is undergoing strong starburst activity. We detect Lya 
emission across ~20 arc sec. Two main components 
are revealed by our spectra, coincident with (LI and 
L2). Lya is relatively broad in LI (FWHM-1800 km 
s- 1 ) and narrower in L2 (FWHM-300-700 km s" 1 ). A 
high velocity region (FWHM~T500 km _1 ) is detected 
at the border of L2. A spectrally unresolved region lies 
at ^12 arc sec from LI. 

The four objects present complex kinematics and show 
that there is a large variety of kinematic behaviour in high 
redshift active galaxies. All objects present certain com- 
mon characteristics: 

- High velocities (FWHM>1000 km s _1 ) in the extended 
gas 

— Velocity shift of the Lya emission across the nebula 
varying over a range <700 km s , although the ve- 
locity curves are rather different from object to ob- 
ject. Similar values are observed in low redshift radio 
galaxies (Tadhunter et al. 1989) 

— Presence of at least two different kinematic compo- 
nents which seem to be spatially distinct. Such com- 
ponents look like individual clumps in the case of 
MRC2025-218, MRC2104-242 and SMM02399-0136. 
Diffuse and fainter line emission is present in the spec- 
tra of SMM02399-0136 and MRC1558-003. Narrow 
band Lya narrow band images show also diffuse Lya 
emission in MRC2025-218 and MRC2104-242 in ad- 
dition to the brightest components (McCarthy et al. 
1990, Pentericci et al. 1998). 

5. Discussion 

High velocities have been observed in the EELR of many 
HZRG (McCarthy et al. 1996). The alignment between 
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Fig. 1. MRC2104-242. Upper panel: 2-D spectrum of the 
Lya spectral region with the dispersion in A running ver- 
tically. Middle: Spatial variation of the FWHM. Lower 
panel: Velocity shift with respect to the emission in the 
center of symmetry of the Lya distribution. The spatial 
shift has been calculated with respect to the same posi- 
tion. The horizontal lines indicate the pixels added along 
the spatial direction. The three panels have the same spa- 
tial scale and are aligned so that vertical lines connect the 
same spatial positions in the slit. The errorbars (vertical 
lines) include the errors of the fit, as well as the error in 
the measurement of the instrumental profile. 



the radio and optical structures (McCarthy et al. 1987, 
Chambers et al. 1987) and the anticorrelation between 
the size of the radio source and the velocity dispersion 
found for HZRG (van Ojik 1995) suggest that the jet is 



Fig. 2. MRC2025-218. Panels and symbols as in Fig. 1. 
The spatial and velocity shifts have been calculated with 
respect to the centroid of continuum emission. 



interacting with the ambient gas. Studies of radio galaxies 
at intermediate redshift with clear signs of such interac- 
tions as well as hydrodynamical simulations show that 
this process can produce large FWHM (>1000 km s _1 ) 
(e.g. Villar-Martm et al 1999, Clark et al. 1997). Some 
HZRG show clear evidence for jet-cloud interactions (e.g. 
van Ojik et al. 1996) and this process is surely having an 
effect in some high redshift radio galaxies. It could be also 
the case of MRC2025-214 and MRC2104-242, which show 
the alignment effect (McCarthy et al. 1990) and large line 
widths inside the radio structures. However, we have also 
measured high velocities in the EELR of MRC1558-003 
beyond the radio structures (Fig. 3) and the extended gas 
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Fig. 3. MRC1558-003. Panels and symbols as in Fig. 1. 
The spatial and velocity shifts have been calculated with 
respect to the nuclear emission. The high velocity compo- 
nent 10 arc sec from the x =0 position is located beyond 
the radio structures. 

of the galaxy-galaxy interacting system SMM02399-236. 
Bremer et al. (1992) reported the detection of extended 
Lya emission in the radio quiet quasar 0055-264 (z =3.66) 
with FWHM^IOOO km s _1 . Jet cloud interactions cannot 
explain the extreme kinematics in these objects. Another 
accelerating mechanism is at work which could also play 
a role in many other HZRG. 

IV98 have proposed that SMM02399-0136 is a system 
in which two companions are interacting: LI (that con- 
tains an active nucleus) and L2. The interaction has in- 
duced starburst activity responsible for the submm and 
weak radio emissions. Another possibility is that the ra- 



Fig. 4. SMM02399-0136. Panels and symbols as in Fig. 1. 
The two components of the system (IV98) have been in- 
dicated (LI and L2). The spatial and velocity shifts have 
been calculated with respect to the centroid of continuum 
emission. Note the presence of large velocities in the ex- 
tended gas of this very weak radio source. 



dio emission is due to a frustrated radio jet, since the 
steep radio spectrum is consistent with an AGN origin 
(IV98). The radio emission is extended along PA71, close 
to PA88.6 (L1-L2 direction) and the projected size (~7.9 
arc sec) is larger than the L1-L2 extension (~4 arc sec). 
Therefore, the high velocity gas detected beyond L2 is 
probably inside the radio structures. In this case, the in- 
teraction between the radio jet and the ambient gas could 
be responsible for the large FWHM values measured be- 
yond L2. 
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However, if SMM02399-0136 is a system of two interac- 
tion companions with strong starburst activity [as is the 
case for many /^Jy radio sources (Lowenthal 1997)], jet- 
cloud interactions cannot explain the gas kinematics in 
SMM02399-0136. Q 

We discuss here several possible mechanisms that 
could explain the extreme kinematics in some HZRG. 

— Broad scattered lines. 

Many distant radio galaxies (z >2) show polarized con- 
tinuum with the electric vector perpendicular to the 
axis of the optical (UV rest frame) structures (Cimatti 
et al. 1998, Fosbury et al. 1998b,1999). This is consis- 
tent with an scenario in which powerful radio gala- 
xies contain QSO nuclei whose FUV emission we see 
scattered by extended dust structures. The emission 
from the broad line region should also be scattered and 
therefore broad lines could be detected within the ex- 
tended gas. Scattering preserves the equivalent width 
of BLR lines against the nuclear continuum and if this 
mechanism dominated, we should measure similar EW 
in the extended gas. We have measured a lower limit 
for the EW of the CIVA1550 line (not affected by neu- 
tral hydrogen absorption). We obtain EW(CIV)>100, 
which is quite large compared with typical values mea- 
sured in quasars at high rcdshift (Corbin & Francis 
1994). This suggests that the line is dominated by di- 
rect light, rather than scattered. 

This mechanism is not likely to play a role in ra- 
dio galaxies in general. In effect, NIR spectroscopy 
of distant narrow line radio galaxies (z =2.2-2.6) re- 
veals FWHM>1000 km s _1 for both permitted and 
forbidden lines (Evans 1998). Eventhough the spectra 
of HZRG are customarily integrated along the spa- 
tial dimension, the extended emission frequently domi- 
nates the integrated spectra (e.g. McCarthy et al. 1996, 
Stockton et al. 1996) and therefore the observed large 
velocities could well originate from the extended re- 
gions themselves. At z ^1 the [Oil] emission has simi- 
larly revealed high velocities (FWHM>1000 km s _1 ) 
within the extended gas regions of some radio galaxies 
(McCarthy et al. 1996). 

— In jail. 

Heckman et al. (1991) made a spectroscopic study of 
the extended gas in a sample of 5 high redshift radio- 
loud quasars (z ~2-3) (see also Lehnert & Becker 
1998). They found that the kinematical properties are 

1 SMM02399-0136 is gravitationally lensed by a foreground 
cluster. The appropriate geometry with respect to the lensing 
source could produce a greater distortion of component LI due 
to its compact morphology. LI emission might extend beyond 
L2 with the result that the high velocities measured in LI could 
be contaminating measurements of the extended gas. We would 
then expect the same effect for NV and CIV, which are strongly 
nucleated and have similar fluxes in LI than Lya (see Fig. 3 
in IV98). However, these lines are detected only in LI. 



very similar to the properties of the EELR of HZRG: 
namely velocity dispersions across the nebulae consis- 
tent with gravitational motions (<500 km s" 1 ) but 
with large FWHM~1000-1500 km s" 1 . They propose 
a scenario in which gravitation is at the origin of these 
extreme motions: i.e. gas freely falling from a large dis- 
tance into the galaxy. An infall process with such cha- 
racteristics could happen during the process of galaxy 
formation, if the radio galaxy lies at the bottom of a 
deep potential well (like a dense cluster). 

— A group of Ly break galaxies around the radio galaxy. 
Recent deep HST images of radio galaxies at z >2 show 
clumpy and irregular morphologies, consisting of a 
bright component and a number of small components, 
which is suggestive of a merging system (Pentericci 
et al. 1998). Pentericci et al. find that those clumps 
have similar characteristics to Ly break galaxies and 
suggest that the host galaxy of the radio source had 
itself formed through the merging of such smaller units 
(see also van Breugel et al. 1998). 
MRC2025-218, MRC2104-242 and SMM02399-0136 
have similar morphologies to the one described above 
(see optical images in Pentericci et al. 1998 and IV98). 
The presence of several components is also revealed by 
our spectra. We also detect diffuse emission between 
(and sometimes beyond) the main clumps which could 
have the same nature as the diffuse and asymmetric ha- 
los found around compact clumps in Ly break galaxies 
(Steidel et al. 1996). 

An important difference between Ly break galaxies and 
the clumps in HZRG is the large velocities measured in 
each clump (FWHM>1000 km s _1 ), much larger than 
the values observed in Ly break galaxies (FWHM<200 
km s _1 ) (e.g. Pettini et al. 1998, Prochaska & Wolfe 
1997). However, spectroscopy shows that velocity 
shifts of >1000 km s _1 between absorbing and emit- 
ting gas are common in Ly break galaxies. This su- 
ggests the presence of large scale outflows of hundreds 
km s _1 (e.g. Pettini et al. 1998, Franx et al. 1997), 
which could be responsible for large FWHM values if 
all the gas was ionized. On the other hand no clear link 
has been established between the absorbing and the 
ionized gas and, therefore, these velocity differences 
might simply occur between unrelated intervening ob- 
jects. 

— Bipolar outflows. 

Chambers (1998) has recently proposed that bipolar 
outflows can be responsible for the high velocities, mor- 
phologies and polarimetric properties of high redshift 
radio galaxies. The EELR of HZRG consists in this 
model of an expanding bipolar dust shell which scat- 
ters light from a quasar core and has an evacuated 
interior. 

Bipolar outflows can be generated by the superwind 
associated with a starburst in a circumnuclear mole- 
cular disk. Evidence for such superwinds has been 
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found in Far Infrared Galaxies (FIRGs) (Heckman et 
al. 1990, H90 hereafter). These galaxies show emission 
lines with FWHM of several hundreds km s _1 and 
shifted by —1000 km s _1 in some objects. Appropriate 
superwind models predict outflow velocities of several 
hundreds km s . 

We have calculated some basic parameters characte- 
rizing a superwind which could explain the kinematic 
properties of the EELR in HZRG (we want to en- 
phasize the approximate nature of this calculations). 
We have assumed a typical radius r ne b for the ioni- 
zed nebula of 20 kpc, a density no in the ambient 
(undisturbed) medium of 10 cm" 3 (McCarthy 1993). 
Two emission lines components of several hundred 
km s _1 (unresolved at our spectral resolution) and 
shifted by 1000 km s _1 will produce a broad profile 
with FWHM— 1000 km s -1 , the values observed in our 
objects. Therefore, we can assume, as concluded for 
nearby FIRGs, an expanding velocity of the nebula 
v ne i, of several hundreds km s _1 (500 km s -1 ). The 
dynamical time scale td yn for the nebula is given by 
equation [7] in H90. We obtain td yn ~40 Myr, which is 
in the range calculated for some FIRGs. According to 
the predictions of the superwind models, this is compa- 
rable to the age of the starburst. For comparison, Dey 
et al. (1997) derived an upper limit for the young ste- 
llar population in the radio galaxy 4C41.17 (2=3.80) 
of 600 Myr for a continuous star formation scenario 
and ~ 16 Myr for an instantaneous starburst model. 
On the other hand, if the total mass of the ionized gas 
(severalxl0 8 -10 9 M Q , van Ojik 1995, McCarthy 1993) 
has been ejected in the outflow and td yn ~40 Myr, this 
implies that the mass injection rate is 7< ^ <25 M Q 
yr~\ which is also consistent with superwind models 
for nearby FIRGs. 

We have also calculated the rate of injection of kinetic 
energy in the wind given by eq. [8] in H90. We 
obtain — 3xl0 47 erg s^ 1 which is one order of mag- 
nitude higher than predicted for high power FIRGs. 
This is what we expect taking into account cosmologi- 
cal evolution of the wind rate (H90) (the mass injection 
rate should increase by the same factor, though). ^ 
is large compared to the integrated luminosity of the 
lines (Lycv can be as luminous as — 10 44 erg s _1 ).With 
a small conversion efficiency of kinetic energy of the 
outflow in emission line luminosity, the outflow could 
power the nebula. 

Therefore, a superwind with properties similar to those 
predicted for the high power FIRGs at low redshift, 
could explain the kinematic properties of the extended 
gas in the HZRG. 

The possibility of an outflowing wind generated in the 
nuclear region of a powerful radio galaxy is suggested 
by the nearby radio galaxy Cyg A: a polarization map 
(Tadhunter et al. 1990, Ogle et al. 1997) shows a bi- 
conical structure suggestive of a dusty reflection nebula 



which scatters the light of the hidden active nucleus. 
The outflow is suggested by the shift between the na- 
rrow emission lines detected in direct light and the 
polarized narrow lines (Ogle et al. 1997). Evidence for 
a circumnuclear starburst ring is presented in Fosbury 
et al. 1998a. The velocity shift however, is much lower 
(-100-200 km s" 1 ) than the FWHM values observed 
at high z. 

6. Summary and conclusion 

We have studied the UV spectra of 3 distant powerful ra- 
dio galaxies and the hyperluminous SMM02399-136 sys- 
tem with the goal of understanding the mechanism re- 
sponsible for the high velocities observed in the extended 
gas of HZRG. Large velocities are found in the extended 
gas of all the objects (FWHM>1000 km s" 1 ). 

Interactions between the radio jet and the ambient gas 
certainly play a role in some radio galaxies. However, we 
measure high velocities in regions where such type of in- 
teractions are not taking place and therefore other mech- 
anisms must be at work. 

Possible explanations for such extreme motions are: 1) 
infall of material from large distances (gravitational ori- 
gin). This mechanism could be important in the process 
of galaxy formation. 2) A group of Ly break galaxies in 
the neighbourhood of the radio galaxy. Large scale out- 
flows in the individual components are required. 3) Bipolar 
outflows produced by superwinds, as observed in nearby 
FIRGs. 
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